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Abstract. We analyse a feature that is detectable at the differential phases of high- 
resolution spectro-interferometry of Be stars. The origin of this feature (dubbed CQE- 
PS, Central Quasi Emission Phase Signature) lies in the differential absorption of pho- 
tospheric radiation by the circumstellar disk, which is spectroscopically detected as a 
CQE line profile in shell stars. This phenomenon has great diagnostic potential for Be 
star disks, revealing properties of these disks on the scale of a few stellar radii. 



1. Introduction 

Classical Be stars are well known to be characterized by having gaseou s circumstellar 



disks that are fed from mass lost from their rapidly rotating central stars (Por ter & Rivinius 



2003). In the past two decades, optical and infrared interferometry played a key role in 



the study of these astrop hysical systems. For instance, the pioneering measurements of 



Ouirrenba ch et all ( 19971) provided strong evidence that the circumstellar matter around 



Be stars is distributed in a geometrically thin structure (i.e., a disk). More recently, de- 
tailed measurements allowed to demons t rate that Be disks rotate i n a Keplerian fashion 
dMeilland et al.ll2007l : iKraus et al.ll20ll IWheelwright et al.ll2012l) . 



The strong developments in the observations were accompanied by significant 
progress in our theoretical understanding of these systems. Recently, a co nsensus is 
emerg ing that the viscous decretion disk model (VDDM), first suggested by iLee etail 
(1199 ll) and further developed by Bjorkman (1997), and Okazaki (2001), among oth- 
ers, is the best cand idate to explain the observed properties of Be disks dCarcionll201ll : 
iMcGill et al.ll20TTh . This model uses the angular momentum transport by turbulent vis- 
cosity to lift material that was ejected by the star into higher orbits, thereby causing the 
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disk to grow in size. A theoretical result of this model is that material is in Keplerian 
rotation throughout the disk, in agreement with the recent interferometric results. 

Before the spectro-interferometric era, one of the first strong indications that Be 
disks are Kepler i an ca me f rom high-reso l ution spectroscopy (R > 30000) of shell 



Hanuschik ( 1995) and Rivi nius et al] ( 1999) studied the so-called Central Quasi- 



stars. 

Emission (hereafter, CQE) peaks, where the disk, under certain circumstances, causes 
a cusp in the deepest region of the line profile of shell stars. The existence of the CQE 
implies slow radial motions of the gas, which means the disk is supported by rotation. 

In this work, we show that the same mechanism that produces the CQE observed 
in shell star line profiles can cause important changes in the intensity map of the Be star 
plus disk system, with observable effects on the interferometric quantities. We analyze 
and quantify the diagnostic potential of this effect, which we dub CQE Phase Signature 
(hereafter, CQE-PS). Interestingly, similarly to the CQE detection, the CQE-PS can 
only be studied by high-resolution observations (R > 12000), such as can be obtained 
by AMBER/VLTI in the near infrared. 



2. Model Description 

We will see in Sect. [3] that the mechanism behind the CQE-PS is the differential ab- 
sorption of photospheric radiation by the circumstellar disk. This mechanism can be 
directly traced by the interferometric differential phase, and must be separated from 
other mechanisms that affect this quantity. 

To study, characterize and illustrate the CQE-PS, we adopt a realistic model for 
the Be + disk system. For the central star we adopt a rotationally deformed and gravity 
darkened star whose parameters are typical of a BlVe star (Table [D. In this study, the 
only free parameter of the star is its rotation rate (Q/Q cr i t ). For the circumstellar disk, 
we adopt the VDDM. For the case of isothermal viscous diffusion in the steady-state 
regime, a state reached after a sufficiently long and stable decretion period (Haubois, 
these proc eedings), the disk sur face density assumes a particularly simple form E(r) oc 
r~ 2 (iBiorkman & Carcionll2005h . 



The interferometric quantities depend on the spatial resolution with which the ob- 
ject is seen. To explore different configurations we define the quantity v t, s , the ratio 
between the distance to the star (i.e., its angular size on sky) and the baseline length of 
the interferometer. The units used are mpc -1 . 



2.1. Interferometric Phases in the Marginally-Resolved Case 

We see below that the CQE-PS at interferometric differential phases (hereafter DP) 
happens under quite specific circumstances, and its diagnostic potential is more relevant 
for some parameters than for others. To quantify this, we first study a reference model 
for which the CQE-PS is not present, i.e., a model for which the canonical signatures 
for interferometric quantities apply. The adopted model parameters for the reference 

case are listed in the thi rd column of Table [J 

As pointed out by iDomiciano de Souza et al.1 (120041) . an image i(r, A) of an astro- 
nomical target is defined as the convolution of the intensity map of the object, I(r, A), 
with the Point Spread Function (PSF) due to the instrument and the terrestrial atmo- 
sphere, A(r, A), i.e., i(r, A) = I(r, A) * A(r, A), where r is the position on the sky vector. 
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Table 1. Reference Be model parameters. 



Parameter 


Symbol 


Ker. ^ase 


Spectral type 


- 


B1V 


Polar radius 


^pole 


4.9 R @ 


Pole temperature 


^pole 


27440 K 


T iimino*;itv 

i—i LiiiiiiivjiJiL y 


i* 


1 01 60 Jr. 
J. \J J. \J\J -I—/Q 


Critical velocity 


v crit 


534.4 km/s 


Rotation rate 


0/Q cr i t 


0.8 


Oblateness 


^eq/^pole 


1.14 


Gravity darkening 


^pole / ^eq 


1.16 


Disk radius 


Rdisk 


10 R« 


Disk density 


n 


10 13 g cm- 3 


Incl. angle 


i 


45° 


Spectral res. 


R 


12000 


Baseline/distance 




0.5 m pc _1 


Position angle 


PA 


0° 



When the data is obtained by an interferometer, the vector u represents the spatial 
frequency associated with r. In long baseline interferometry, u is the vector B pro j//l. 
At high spectral resolution the interferometric phase signal <p can be obtained as an 
expansion of /(u, A), the Fourier transform of the brightness distribution. 

If the target is not r esolved u r« 1, i.e., very small values of v b s > the differential 
phase can be written as (|Domiciano de Souza et al. 2004, Eq. 4) 



0(r, A, A r ) = -2nu ■ [e(A) - e(A r )] ; (1) 

the e(A) and e(A r ) vectors are called photocenters respectively for A and A r , where A r is 
a reference channel observed simultaneously. The usual procedure is to take A r on the 
adjacent continuum of the spectral line, with e(A r ) = 0. 

For the case of rotating disks, it is usually assumed that the DP will have a simple 
S-shaped profile (ISteelll996l) . This can be easily understood from Fig. [T] which shows 
the Br-y line profile, the corresponding interferometric quantities and model images at 
different wavelengths for the reference model. Shown are results for i = 45°. At con- 
tinuum wavelengths (d, in Fig. [TJ) the images show the star surrounded by a centrosym- 
metric continuum emission. We define the DP to be zero at those wavelengths. As we 
move from the continuum towards line center (d to a) the line flux initially increases 
as a result of progressively larger emission lobes and then decreases as the emission 
area decreases towards line center (a). This creates the familiar W-shaped pattern in the 
visibilities and the S-shape pattern for the DPs. 

2.2. Interferometric Phases in the Resolved Case 

The association between the phase signal and the image's photocenter position is only 
possible if the marginally -resolved condition is satisfied, since it permits the truncation 
of the Fourier expansion (iJankov et al.ll200lh . This is particularly important for Be stars 
since their on-sky size is wavelength dependent, mainly across emission lines where the 
size of the emitting regions varies with wavelength. 
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Figure 1 . Top 3 panels: synthetic line flux, interferometric visibilities and phases 
across the Br-y line for the reference model. Calculations were made for v bs = 
0.5 mpc , and baseline orientation parallel to the disk equator. Bottom panels: The 
Be model images for the reference model associated with its high-resolution inter- 
ferometric measurements. The disk inclination angle is i = 45°. 
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Figure 2. Differential phase signal for the Be reference model. The reference 
model is the black (full) line, observed with v b s = 0.5 mpc -1 (astrometric regime). 
In the red (dashed) line, v t, s = 2.5 mpc -1 (complex phase behaviour). 



Outside the astrometric regime (i.e., u • r ~ 1 or large values of v hsX Eq. © is 
no longer valid. A complex phase behavior is then expected since higher order terms 
of the Fourier transform of the image are relevant, and it does not correspond anymore 
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Figure 3. The squared visibility and differential phases as function of inter- 
ferometric baselines for the Be star reference model at radial velocities of v = 
-lOOkms' 1 and v = -152kms _1 . 



directly to the photocenter. Fig.|2]compares the DP signal of our reference model in the 
astrometric regime and outside it. The reason the curves are so different is explained 
in Fig. [3] that plot the interferometric quantities for two wavelengths accross Br-y as a 
function of baseline length. The complex phase behavior, even with change of sign, 
occurs even before the fully resolved regime, i.e., when the squared visibility goes to 
zero. 

Based on our simulations, the astrometric regime is approximately valid up to 
squared visibilities of « 0.8, which corresponds to v b s < 1.5 mpc" 1 for a typical Be 
star. The exact value of these quantities is dependent of the target brightness distribu- 
tio n. One example of the complex behavior outside the astrometric regime is reported 
bv lKraus et ail (12012b . who models AMBER/VLTI data for the Be star /? CMi. In this 
case, a phase reversal is registered at the DP signal of the longest baseline available for 
this nearby target (squared visibility of « 0.7, v b s * 2.2). 

Outside the astrometric regime a phase reversal is an example of the phenomenol- 
ogy that emerges entirely from long baseline interferometry, where the phase signal no 
longer corresponds to displacements of the photocenter. In the text that follows all inter- 
ferometric computations are in the marginally resolved regime to adequately separate 
the CQE-PS phenomenon from the just described effects. 



3. The CQE Phase Signature 

For the reference model described in Sect. 12.11 line emission from the rotating disk is 
the most important factor controlling the interferometric quantities. However, in some 
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circumstances line absorption of photospheric light by the disk can also have a strong 
impact on these observables. 

The spectroscopic phenomenon of the CQE was first studied bv lHanuschikl(ll995h 
for the so-called shell absorption lines in nearly edge-on Be shell stars. These consist of 
narrow absorption lines well below the stellar continuum that show a cen tral weak cusp, 
known as Central Quasi-Emission peaks (CQEs). As pointed out by iRivinius et all 
(1999), these peculiar profiles result from the minimum at zero radial velocity in the 
absorption of the photospheric flux by the rotating circumstellar gas. 

CQEs are, in spite of their name, not related to any emission process but are a pure 
absorption phenomenon. So, when observing the star plus disk system the latter acts as 
an absorber of the stellar photospheric light, as seen in the CQE spectroscopy profile. 

The differential absorption across the spectral line that causes the CQE also has 
strong effects on interferometric differential phases, as shown in Fig. |4j where we plot 
results for the reference case seen at i = 90° (edge-on). The projected absorption in 
the stellar disk generates a considerable decrease of the stellar flux, which in turn will 
affect the photocenter of the system, with a corresponding signal in the DP. 

Consider the system seen at continuum wavelengths (Fig. |4ji). The dark lane 
across the stellar disk is caused by continuum (bound-bound and bound-free) absorp- 
tion. As we consider images at different spectral channels across the Br-y line, we 
also see line (bound-bound) absorption of photospheric light and emission by the disk 
(Figs. to|4fc). This line absorption has a differential aspect because it depends on the 
line of sight velocity of the absorbing material. 

At high velocities, only a small area near the edge of the disk will block the 
starlight (Fig. |4fc). For lower velocities, an increasing disk portion will block the 
starlight, creating a maximum absorption around 50 km s _1 (Fig. §J)). Finally, as one 
considers wavelengths closer to the line center, the disk absorption decreases due the 
fact that there is less disk material with the appropriates line-of-sight velocity. 

The CQE effects are concentrated very near the central wavelength, since there the 
opacities along the line-of-sight are largest. Fig. [5] shows the projected velocity interval 
of the disk as a function of the stellar impact parameter. The low-velocity components 
are present for all impact parameters, and that is why the CQE-PS affects more the DP 
close to line center. 

Contrary to disk emission that shifts the photocenter position towards the direction 
of the emitting lobe, line absorption of photospheric light shifts it on the opposite di- 
rection. So, when evaluating the differential phases, we are measuring the competition 
between the photocenter shift due to disk emission vs. the shift due to disk absorption. 
If absorption is strong enough, or emission is weak (for instance, as expected for a small 
disk density) one may observe a reversal in the DP. 



4. Conclusions 

The new interferometric observations of Be stars at high spectral resolution lead to new 
measurements of differential phases. Deviations from the canonical appearance may 
occur in two different circumstances: (i) by purely interferometric properties at obser- 
vations outside the marginally resolved regime, where the differential phases no longer 
correspond to the photocenter of the system, and (ii) by physical process accessible at 
high-resolution affecting the interferometric signal. In particular, in this work we dis- 
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Figure 4. Top 3 panels: synthetic line flux, interferometric visibilities and phases 
across the Br-y line for the reference model. Calculations were made for v bs = 
0.5 mpc , and baseline orientation parallel to the disk equator. Bottom panels: The 
Be model images for the reference model associated with its high-resolution inter- 
ferometric measurements. The disk inclination angle is i = 90°. 



cuss the role of the disk differential absorption as an important factor in shaping the 
differential phases of Be stars. 

The solution of the radiative transfer problem of stellar photospheric flux in inter- 
action with a typical rotating Be disk can cause a measurable photocenter displacement. 
Since this is an effect of differential absorption by the disk, it will depend on various 
system parameters and can provide valuable information about its configuration. The 
diagnostic potential of the CQE-PS in determining important disk parameters, such as 
density and velocity field, will studied in a future publication (Faes et al., in prep.). 
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Figure 5. Projected velocity range of a rotating disk as a function of impact pa- 
rameter. The solid line shows the average geometric velocity component for each 
impact parameter for a Keplerian disk of = 10/?+. Disk inclination angle is 
i = 90°. 
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Discussion 

Millour: Your "CQE effect" is a physical effect which should be seen at all baselines 
lengths, including spectro-astrometric measurements on Be stars. On the other hand, a 
resolution effect (spatially speaking) would produce a "phase reversal" which depends 
on the baseline length, i.e. it would be seen only on longest baselines and not on 
spectro-astrometry data. My experience would make me point to that hypothesis. What 
is yours? 

Faes: Indeed long baselines lead to complex phase signals - including a reversal. We 
were aware of this in our simulations. The phase reversal presented is assigned to 
photocenter changes from a physical cause, namely stellar light absorption, that occurs 
for specific circumstellar disk configurations. 



